Nonequilibrium roughening transition by two-species particles 
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We introduce an interface growth model exhibiting a nonequilibrium roughening transition 
(NRT). In the model, particles consist of two species, and deposit or evaporate on one dimen- 
sional substrate according to a given dynamic rule. When the dynamics is limited to occur on 
monolayer, this model has two absorbing states, belonging to the directed Ising (DI) universality 
class. At criticality, the density of vacant sites at the bottom layer in the growth model decays 
faster than the DI behavior, however, the dynamic exponent is close to the DI value, suggesting 
that the dynamics is related to the DI universality class. We also consider an asymmetric version 
of the growth dynamics, which is according to the directed percolation behavior. 

PACS numbers:05.70Fh,05.70Jk,05.70Ln 



Phase transitions in nonequilibrium systems are dis- 
tinct in their physical properties from ones in equilibrium 
systems, in general Q. For example, roughening transi- 
tion (RT) from a smooth phase to a rough phase in one 
dimension does not occur in thermal equilibrium systems 
0. In nonequilibrium systems, however, there are a few 
examples of exhibiting RT even in one dimension, such 
as the polynuclear growth models ||, the deposition- 
evaporation model with no evaporation on terrace [Q, 
and the fungal growth model g. The common feature 
of the models is that some aspects of the RT are re- 
lated to directed percolation (DP) in 1+1 dimensions Q. 
On the other hand, DP behavior occurs in many other 
nonequilibrium systems exhibiting phase transition from 
active to inactive state 0. Typical examples would be 
the monomer-dimer model for the catalytic oxidation 
of CO g, the contact process (9), the surface depin- 
ning models [ p)| |, and the branch-annihilation random 
walks with odd numbers of offspring fTl|| . Recently a 
lot of efforts have been made to search for other sys- 
tems exhibiting non-DP behavior. As a result, a few 
examples have been found such as the probabilistic cel- 
lular automata model jjj], certain kinetic Ising model 
p3| , the interacting monomer-dimer model |l4| , modified 
Domany-Kinzel model p5[ , and the branch-annihilation 
random walks with even numbers of offspring (BAWe) 
JyJ . All these models except for the BAWe have two 
equivalent absorbing states indicating the importance 
of symmetry in absorbing states to classify universal- 
ity classes. Analogous to equilibrium spin models, the 
non-DP class with two equivalent absorbing states is re- 
ferred to as the directed Ising (DI) universality class jl6) . 

Recently the stochastic models [3-5] exhibiting 
nonequilibrium RT (NRT) were introduced, all of which 
do not possess any absorbing state. When absorbing 
state is not present, NRT was understood in terms of 
the spontaneous symmetry breaking of non-conserved 
order parameter in Refs.[4-5]. In this paper, we consider 
the classification of universality classes in NRT having 
no absorbing state. Analogous to phase transition with 
absorbing state, we think that the symmetry lying in 



dynamic rule can replace the role of symmetry in ab- 
sorbing states, and characterizes universality classes. In 
order to confirm this idea, we study a simple deposition- 
evaporation model exhibiting NRT, in which particles 
consist of two species, and its dynamic rule is symmetric 
with respect to particle species. We also consider an 
asymmetric version by breaking the symmetry in the 
dynamic rule. Performing numerical simulations, it is 
found that the symmetric and the asymmetric cases be- 
have completely distinctively, confirming that NRT is 
characterized by the symmetry of dynamic rule. Fur- 
thermore, we investigate the classification of universality 
classes for the symmetric and the asymmetric case, re- 
spectively. For the symmetric case, the density of vacant 
sites at the bottom layer at the threshold of NRT, de- 
cays faster than the DI behavior, however, the dynamic 
exponent appears close to the DI value, suggesting that 
the growth dynamics is still related to the DI universal- 
ity class. For the asymmetric case, the density behaves 
well according to the DP class. We will discuss physical 
reason of the different behaviors later. 

The stochastic model we introduce is defined as fol- 
lows. We consider an interface dynamics with deposition 
and evaporation processes of particles on one dimen- 
sional substrate with the periodic boundary condition. 
Particles consist of two species, black and gray colored. 
The dynamics starts from vacuum state, where no par- 
ticle is in the system. First, a site is selected at random. 
Next, either deposition or evaporation of a particle is 
attempted at the selected site with probability p (p/2 for 
black and p/2 for gray particle) and 1 — p, respectively, 
and the attempt is realized under the two conditions be- 
low. First, a restricted solid-on-solid (RSOS) condition 
is imposed such that height difference between nearest 
neighboring columns does not exceed one. Second, we 
consider a ferromagnetic interaction between two nearest 
neighboring particles in the same layer. The interaction 
is attractive (repulsive) when they are the same (differ- 
ent) colored. Then a particle with a certain color cannot 
deposit (evaporate) when both of neighboring particles 
in the left and the right sides have different (the same) 



color from its own one as shown in Fig.l. Particle can 
deposit or evaporate when the colors of two neighboring 
particles on each sides are alternative, or one (or both) of 
the neighboring sites is (are) vacant. Symmetric dynamic 
rule with respect to particle-color generates symmetry in 
surface configurations; every surface configuration has its 
corresponding configuration with replacing one species 
by the other species as shown in Fig.l. 

In an asymmetric version, one-species particles, say 
black colored, can distinguish the same species from the 
other species, however, gray particles cannot see colors 
and regard black particles as the same species. Then 
a gray particle can deposit (cannot evaporate) even at 
the site where both of two neighboring particles in each 
sides are black. Then the symmetry in the dynamic rule 
breaks up, and the gray flourish more than the black. As 
an extreme case of the asymmetric version, when par- 
ticles are of single species, the model is reduced to the 
one by Alon et al 0, where particles cannot evaporate 
at any site in terrace, and only can do at the edge of the 
terrace. The model defined so far is called the growth 
model hereafter to distinguish itself from the model be- 
low, called the monolayer model, confined in monolayer. 

When the dynamics is restricted on monolayer, so that 
particle cannot deposit on top of particle, the model ex- 
hibits a phase transition from active state to absorbing 
state. In this case, vacant site means active site, and 
site occupied by black or gray particle means inactive 
site. There are two types of inactive sites according to 
colors. Absorbing state is formed when the entire system 
is filled with a single species of particles. Then there 
exist two absorbing states which are equivalent as long 
as the dynamic rule is symmetric with respect to particle 
species. Thus, it is expected that the phase transition 
of the monolayer model is in the DI universality class. 
Our monolayer model is similar to the generalized contact 
process proposed by Hinrichsen [fj"5| , however, our growth 
model is much simpler because it includes one control- 
parameter rather than two. Also, our model is easily 
generalized into the cases of higher symmetry involving 
g-species particles or of higher dimensions. Moreover, 
our model might be relevant to describe the deposition- 
evaporation process of Ising particles on surface. 

For the growth model, when p is small, a smooth phase 
is maintained. In this phase, particles form small-sized 
islands which disappear after their short lifetime. Thus 
the growth velocity is zero in the thermodynamic limit. 
As p increases, deposition increases and typical islands 
grow, until, above a critical value p Cl islands merge and 
fill new layers completely, giving the interface a finite 
growth velocity. Accordingly, surface exhibits NRT from 
a smooth phase to a rough phase across p c . We measure 
the density of vacant sites p g (p, t) at the bottom layer 
averaged over all runs, where the subscript g means the 
growth model. p g (p, t) is saturated at a finite value for 



p < p c and decreases to zero exponentially for p > p c in 
the long time limit as shown in Fig. 2 (a). At criticality, 
p g (p c ,t) scales algebraically as 

p g ( Pc ,t)~t-^. (1) 

We performed Monte Carlo simulations for various sys- 
tem sizes L — 10 ~ 1000, and examined the power- 
law behavior by scanning the probability p. For the 
asymmetric case, p c is estimated as p c « 0.3796, and 
P/v\\ ~ 0.16(1) is measured, which is in good agreement 
with the DP values, f3/v\\{DP) « 0.1595 [§. 

For the symmetric case, p c is estimated as p c « 0.4480. 
P/v\\ ~ 0.58(1) is measured in the long time limit, which 
is much deviated from the DI value, (3/v\\(DI) w 0.27 ~ 
0.29 [12-17]. We think that this discrepancy may come 
from the suppression-effect by particle on upper layer; 
For example, particle A in Fig.l can evaporate in the 
monolayer model because one of neighboring particles 
has different color, however, it cannot evaporate in the 
growth model because of another particle on top of it. 
On the other hand, for the asymmetric case, particle 
A cannot evaporate even in the monolayer model, be- 
cause two particles on each side are regarded as the 
same species. Therefore, the suppression-effect appears 
strongly (weakly) for the symmetric (asymmetric) case. 
When number of particles on upper layers is relatively 
small for p <C p c , or in the short time regime for p = p c , 
the suppression-effect is relatively weak. Thus, p g (p c ,t) 
decays according to the DI behavior in the short time 
regime as indicated by dashed line in Fig. 2 (a). However, 
in the long time limit at p c , evaporation is suppressed by 
particles on upper layer, and p g (p Cl t) decays as ~ <~°- 58 , 
faster than the DI behavior. Note that the value * 0.58 
is almost double to the DI value. 

We consider the probability S(t) that the system con- 
tains at least one vacant site on the bottom layer, and 
other sites are occupied by any species of particles. Note 
that S(t) is different from the survival probability P(t) in 
the monolayer model, in that the probability S(t) (P(i)) 
counts the configurations that the bottom layer is filled 
with particles of any species (a single species). As can 
be seen in Fig. 2(b), there exists a characteristic time r g 
such that for t < r g , S(t) — 1, and for t > r g , S(t) decays 
in the same way as does p g . Thus the density of vacant 

sites p g (p c ,t), averaged over samples with at least one 
vacant site, decays as Eq.(l) until r s , and is finite be- 
yond T g as shown in Fig. 2(c). The steady state value of 
Pg S \p c ) depends on system size L as p g s \q c ) ~ L^ 13 ^^. 
We obtained /3/^j_ ~ 0.95(1), which is almost double to 
the DI value, w 0.5. However, the value of the dynamic 
exponent we measured v\\/v_\_ ks 1.64(2), describing the 

characteristic time via r g ~ L v \^ 1/± , is close to the DI 
value, w 1.66 ~ 1.75. We also examine the steady-state 
density p g (j>, oo) for various e = (p c — p) > 0, which may 
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be written as p g (p, oo) ~ e* 3 . As shown in Fig. 2(d), the 
data do not fit well to a straight line for small e, but 
are likely to approach a line asymptotically with slope 
(3 « 0.88, the DI value, for large e, being far from p c . 
Accordingly, the dynamics of the growth model seems 
to be related to the DI class, but the correlation length 
exponents u\\ and vj_ seems to be halves of the DI values. 

In the rough phase, surface grows with finite velocity, 
which behaves as v ~ (q — q c ) v . The velocity would be 
characterized by the inverse of the characteristic time 
T g to complete one layer, v ~ a/r g , where a is lattice 
constant. If r g is regarded as the characteristic time r m 
for the monolayer model to reach the absorbing state, 
T m ~ (q — ?c) -!/ ") where the symbol of bar means the 
monolayer model, then the velocity exponent y would be 
equal to Du . This assumption is confirmed in the asym- 
metric case by that the value measured, y w 1.67(8), 
is close to the DP value, v\\{DP) sa 1.73. However, for 
the symmetric case, we measured 2y pa 2.50(6), which is 
deviated from v\\{DI) pa 3.17 ~ 3.25 @. This discrep- 
ancy may result from that the suppression-effect appears 
much strongly in rough phase. 

Next, we also consider the surface fluctuation width, 

W 2 (L,t)_ = £M(f) - (£EiM*)) 2 s ^re h t {t) de- 
notes height at site i at time t. The interface width be- 
haves as 

W 2 (Lt)~ / <2C/Z ' fori«L 1 /-, 

(£><) \(-e)*£ a C, fort>£V*, ^ 

where C and z are called the roughness and the dynamic 
exponent, respectively, x is the exponent describing 
the vanishing of the roughness as the transition is ap- 
proached. The exponent \ i s measured as \ ~ 0.34(1) 
(« 0.89(4)) for the symmetric (asymmetric) case as 
shown in Fig.3(b). The asymmetric value is close to 
the value 0.92 obtained from the single-species model jjj. 
The roughness exponent ( is measured to be £ « 0.50(1), 
which is very close to the value of the Edwards- Wilkinson 
(EW) universality class Ml and the Kardar-Parisi-Zhang 
(KPZ) universality class |19| in one dimension. The value 
of ( pa 0.5 is also obtained from the height-height cor- 
relation function, C 2 (r) =< (h(r) — h(0)) 2 >^ r 2C > in 
the long time limit as shown in Fig. 3(c). The growth 
exponent £/z exhibits a broad range of numerical values 
from 1/4 to 1/3 as the probability p increases as shown 
in Fig. 3(d). Since the velocity for p > p c is nonzero, 
the EW universality class for surface growth has to be 
excluded. Nevertheless, since the velocity is extremely 
small near p c , the EW behavior appears numerically for 
finite system, but in the thermodynamic limit, surface 
growth for p > p c is expected to belong to the KPZ class. 
At p Cl the roughness of surface exhibits the marginal be- 
havior, W 2 ~ logt for t <C L x l z and W 2 ~ logi for 
t > L 1/z . 



Next, we perform numerical simulations for the mono- 
layer model. We start with two distinct initial configura- 
tions. First, the initial configuration includes one active 
site, so that all sites are occupied by a single species 
particles except one empty site. Then number of active 
sites increases with increasing time. We measure the sur- 
vival probability P(t) (the probability that the system 
is still active at time t), the density of active sites p m (t) 
averaged over all runs, and the mean-square-distance 
of spreading of the active region R 2 (t) averaged over 
surviving runs. At criticality, p c pa 0.7485, the values 
of these quantities scale algebraically in the long time 
limit as P(t) ~ H, p m (t)_~ f>, and R 2 {t) ~ t* . We 
measured the exponents as S ~ 0.287(1), fj k 0.000, and 
z w 1.141(3), all of which are i n g ood agreement with 
those in the DI universality class U%. Next, we start with 
another initial configuration where all sites are empty. 
This initial configuration is considered for the purpose 
of comparing with the growth model. In this case, num- 
ber of active sites decreases with increasing time. The 
density of active sites behaves as p' m (t) ~ t~ v ' and the 
survival probability does as P'(t) ~ t~ s with the expo- 
nents rj' w 0.279(1) and 5' w 0.000, respectively. Even 
though the two different initial configurations yield the 
different values of the exponents, their sum S + fj seems 
to be the same as 5' + fj', describing the growth of the 
number of kinks over surviving samples in the interacting 
monomer-dimer model |l4| ]. Detailed numerical results 
for the growth model and the monolayer model will be 
published elsewhere p0[ . 



In summary, we have introduced an interface growth 
model with deposition and evaporation of two species 
particles, exhibiting a nonequilibrium roughening tran- 
sition (NRT) from a smooth phase to a rough phase as 
deposition-attempt probability p varies. The dynamic 
rule of the deposition and evaporation process is assigned 
symmetrically and asymmetrically with respect to parti- 
cle species, respectively. For the two cases, NRTs be- 
have distinctively from each other, suggesting that the 
universality class of NRT is characterized by the sym- 
metry lying in dynamic rule. For the symmetric case, 
some features of NRT are related to the DI class, but the 
correlation length exponents Vu and is± in the growth 
models are smaller by halves than those of the mono- 
layer model in the DI class, which may result from the 
strong suppression-effect by particles on upper layer. For 
the asymmetric case, however, NRT are described by the 
DP behavior. 
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FIG. 1. The dynamic rule for deposition (a) and for evap- 
oration (b) for the symmetric case of the growth model. 
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FIG. 2. (a) Double logarithmic plot (DLP) of p g (p,t) ver- 
sus time t for probabilities p = 0.4410 (top), 0.4480(= p c ) and 
0.4600 (bottom). The data are obtained for L — 100, aver- 
aged over more than 500 configurations. The dashed and the 
dotted lines have slope 0.276 and 0.589, respectively, drawn 
for the eye. (b) DLP of S(t) versus time t for the same cases 
as (a), (c) DLP of p\ , (p c ,i) versus t for different system 
sizes L = 50(top),100, 200 and 500. The dotted line has slope 
0.589, drawn for the eye. Inset: DLP of r g versus i at p c . 
(d) DLP of p g (p, oo ) versus e. The dashed line has slope 0.88, 
drawn for the eye. All plots (a-d) are for the symmetric case 
of the growth model. 
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FIG. 3. (a) Plot of the velocity 
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versus — e for the 
symmetric case of the growth model. The dashed line 
v — 0.95 * (— e) 1 ' 26 was obtained by a least-square- fit. (b) 
DLP of W 2 versus — e to measure the exponent \- The dot- 
ted line, a guideline to the eye, has slope 0.34. (c) DLP of the 
height-height correlation function versus distance r at p = 0.9. 
The data are well fit to a straight line with slope 2£ = 1. (d) 
DLP of W 2 versus t at probabilities p — 0.5 and 0.9(top). 
The dotted and dashed lines have slopes 0.49 and 0.62(top) 
to guide to the eye. All plots (a-d) are for the symmetric case 
of the growth model with L = 500, and the data are averaged 
over 500 configurations. 
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